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■ ' The experimental investigation of superfluid turbulence in ^He-B is generally 

not possible with the techniques which have been developed for '^He-II. We 

'^ . describe a new method by which a transient burst of turbulent vortex expan- 

5 i sion can be generated in ^He-B. It is based on the injection of a few vortex 

Qj ' loops into rotating vortex- free flow. The time- dependent evolution of the 

quantized vorticity is then monitored with NMR spectroscopy. Using these 

techniques the transition between regular (i.e. vortex number conserving) 

^ ' and turbulent vortex dynamics can be recorded at T ^ 0.6 Tc and a number 

f— s ■ of other characteristics of turbulence can be followed down to a temperature 

(N: o/r<0.4Tc. 

^ : PACS numbers: 47.37, 67.40, 67.57. 
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INTRODUCTION 



F^ ■ Quantum turbulence, the seemingly chaotic motion of quantized vortex 

P^ \ lines in a disordered network, is a prominent characteristic of superfluid He- 

^ ■ 11, where it has been known to exist for almost half a centuryPUp to recently, 

- - He-II was the only superfluid system in which this phenomenon has been 

investigated. ^He-II is believed to display turbulence in the entire superfluid 
/\ ' temperature range which so far has been probed with measurements. In 

j^ ■ superconductors friction in vortex motion is always high and turbulence is 

not observed: Only in the extreme clean limit one might hope to see any 
signs of other behavior. At present the closest superconducting analogues 
are dynamically driven vortex avalanches."^ 

Superfluid ^He is an intermediate case between these two extremes. The 
dynamics of vortex motion in the ^He superfluids has primarily been studied 
in rotating cryostats. In the past the generally accepted view was that 
the high mutual friction leads to exponentially damped motion and forces 
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vortices to evolve rapidly along well defined trajectories, so that the number 
of vortex lines remains conserved in dynamic processes. A large number of 
measurements over the years has proven this to be the case in both ^He-A 
and ^He-B at T > 0.6 Tc.^ 

Even at low temperatures (~ 0.2 Tc) a duplication of the classic Vinen- 
vibrating-wire measuremenlF' verified that in '^He-B stable vortex motion 
persists for hours in the absence of applied flow. In this experiment super- 
fluid circulation was trapped around a thin wire suspended along the sym- 
metry axis of a cylindrical container while the cryostat was rotating. When 
rotation was stopped, the trapped circulation started to unwind while a vor- 
tex filament, stretched between the wire and the cylinder wall, precessed 
around the wire and spiraled down along the whole length of the 15 mm 
long cylinder. The precessing spiral motion lasted for more than 35 h and 
the precession frequency remained constant with a precision of 0.5%. 

The first indications of rapid non-linear vortex proliferation in ^He-B 
came from measurements of the critical velocity of vortex formation at tem- 
peratures < 0.6 Tc (Ref.ini see p. 140-141). Instead of the usual single- vortex 
formation at a reproducible critical velocity, sudden avalanches of vortices 
were observed in accelerating rotation. These bursts were interpreted to 
mean that at lower temperatures mutual-friction damping had decreased 
sufficiently so that turbulence became possible. 

The first evidence for tangled vorticity was reported from vibrating wire 
measurements in a non-rotating cryostat at temperatures below 0.2Tc.'^A 
vortex network was found to be produced with a vibrating wire resonator 
when it was driven at high level above a critical velocity of vortex formation. 
This was inferred by studying the damping of a second wire, vibrating at low 
level and probing the density of quasiparticle excitations which are Andreev 
retro-refiected by the superflow fields in the vortex tangle. Similar later 
measurements have allowed a determination of the vortex density in the wire- 
generated network. At zero externally applied pressure and at a temperature 
of 0.12 Tc the tangle localized around the vibrating wire turned out to have 
extremely low density, corresponding to an average inter-vortex distance of 
0.2 mm, and to decay away rapidly in about 5 s, when the generator wire 
was switched oS}^ 

There exist no immediately obvious techniques by which turbulence can 
be easily generated and detected in '^He-B. The above two examples suggest 
two ways which are supplementary, since they work in different temperature 
regimes. As outlined below in Sec. EH NMR on a rotating sample can be 
used at temperatures down to about 0.4 Tc. In the zero temperature limit, 
where the density p^ of the normal component approaches zero, moving ob- 
jects can be used to generate vortices and, when driven in a high-quality 
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resonance mode, to detect them. In this report we outUne the techniques 
by which turbulence is generated and detected in the intermediate temper- 
ature regime 0.4 ~ 0.6 Tc, where mutual friction damping and pn are still 
finite. Recently a practical technique was discovered for injecting vortex 
loops into vortex-free flow."^ It then became possible to monitor the nature 
of vortex dynamics as a function of temperature and flow velocity in a con- 
trolled fashion.' ' The transition, which separates regular vortex motion from 
turbulence, could now be mapped. The foundation, on which the interpre- 
tation of these results is built, was laid down by an earlier set of important 
hydrodynamic measurements-^ on mutual friction which provided both the 
dissipative and reactive mutual friction coefficients a(T,P) and a'(T,P) as 
a function of temperature (T) and pressure (P). It turns out that in ^He-B 
the variation of mutual friction with temperature happens to be such that 
the transition falls in the middle of the experimental temperature range, to 
0.5- 0.6 Tc. 

These techniques have been employed in three recent studies investi- 
gating turbulence in '^He-B.mE2I12I xhe firslP^ reports on the flight time of 
the vorticity when it expands along the rotating column. The seconcP^ de- 
scribes how turbulence starts in rotating vortex-free flow in the presence of 
neutron radiation. The thircP*' discusses the mutual-friction dependence of 
the transition from regular to turbulent vortex dynamics, i.e. the transition 
line as a function of T and P at high flow velocities, where the transition 
is velocity independent. These transition properties will also be illuminated 
in Sec. I^of this report, but from a different point of view. NMR measure- 
ments from the transient period, when the turbulent vorticity decays into 
the equilibrium rotating state of rectilinear vortex lines, will be analyzed in 
a forthcoming report. 

According to current view no similar kind of turbulence as in ^He-II 
or ^He-B exists in the anisotropic superfluid ^He-A. Here vortex motion is 
highly damped at all experimentally accessible temperatures and fast vortex 
motion is supported by the appearance of a different structure of vorticity, 
the vortex sheet .-^ 



2. ROLE OF DAMPING IN SUPERFLUID TURBULENCE 

Let us consider a hydrodynamic equation which describes the evolution 
of quantized vorticity in a limit, which becomes useful in the interpretation 
of our measurements. We apply dimensional analysis to this equation in the 
same manner as one characterizes different flow regimes with the Navier- 
Stokes equation in viscous hydrodynamics. This immediately leads to the 
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Fig. 1. Different states of uniformly rotating superfluid. (Left) The vortex- 
free state is the highest energy state and the initial state in these measure- 
ments, with the maximum applied counterflow velocity at given rotation Q: 
V = Vn — Vs = ^f- (Middle) The equilibrium vortex state is the lowest en- 
ergy state. Here an array of rectilinear vortex lines fills essentially all of 
the container and the superfluid component is on an average in solid body 
rotation, i.e. stationary in the rotating frame. The number of vortex lines is 
iVcq ~ vri?^ 20./ K ~ 10'^ if ri ~ 1 rad/s. This is the final state after injection 
of seed vortex loops, if the vortex expansion has been turbulent. (Right) 
A meta-stable vortex cluster may include any number of rectilinear vortex 
lines, < A^ < Agq. Within the central coaxial vortex cluster the areal 
density of the lines 20./ k corresponds to solid body rotation. The counter- 
flow velocity outside the cluster is v = K{Neq — N)/{27rr), where we define 
N = t:E? 20v//v. This is the final state after injection, if loop expansion has 
been regular and the number of individual vortices has been conserved. In 
this case A < 30 after a Kelvin-Helmholtz injection event (Fig. EI). 



notion that a change in the nature of vortex motion is expected as a function 
of mutual friction damping. The argument is simple-minded, but makes the 
suppression of turbulence by mutual friction damping understandable in 
superfluid hydrodynamics. 

The equation of motion for a viscous incompressible fluid is the Navier- 
Stokes equation. 



(9v 
'di 



+ (v • V)v = F/p - VP/p + z^Av , 



(1) 
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where F is an external force per unit volume, P pressure, and i' kinematic 
viscosity. In order to compare flows with different viscosities it is useful to 
write the Navier-Stokes equation in dimensionless units x" = x/L, v^ = v/V 
and t^ = Vt/L, where L and V are the characteristic length and velocity 
scales. The relative importance of the inertial term (v • V)v ~ V'^/L and 
the viscous term uAv ~ uV/Lp' is used to characterize the flow. The ratio 
of these two terms is called the Reynolds number Re = iV'^ / L) / {vV / L?') = 
VLjv. For flows with small Re the dissipative forces are dominant and 
the resulting flow is typically laminar. If Re is large, the inertial forces are 
dominant and the flow is usually turbulent. 

The same type of classification in superfluid hydrodynamics can be 
based on an equation for the "coarse-grained" superfluid velocity Vg. It 
is obtained by averaging over volumes containing many vortex lines. The 
equationP^ is constructed from the Euler equation: 

^ + (vs • V)vs = F/ps - V/x, (2) 

which can also be written as 






,2 



+ c^xvs = F/ps-V(/i + ^j, (3) 



where /x plays the role of the chemical potential and a; = V x Vg is the 
vorticity. 

For such coarse-graining to be meaningful the volume, over which we av- 
erage, should contain many roughly similarly oriented vortex lines. As each 
vortex carries a fixed quantized circulation k = 6.61 • 10~^ cm^/s the number 
of vortex lines in the sample can be estimated from below as [{uj) / k\LP' ~ 
VL/k. We call this combination the "superfluid Reynolds number" Res = 
VL/k and require that Res ^ 1 for the coarse-grained equation Q to be 
applicable. 

The Euler equation basically expresses the conservation of energy and 
Eq. (jSJ thus accounts for the kinetic energy of the global flow. The vortex 
line tension, which is associated with the 1/r local velocity profile around 
the vortex cores, is neglected. The tension can be included by introducing a 
correctiorP^^ to the superfluid velocity dvs = (K/47r) ln(6/a)V x a>, where 
b is roughly the intervortex spacing, a is the vortex core size, and u is the 
unit vector in the direction of u. Dimensional arguments give the estimate 
Vs/Svs ~ V/{k/L) = Res- Thus for values Res ^ 1, the tension can be 
neglected. 

In the classical equation (^ the term F represents some force applied 
externally to the fluid. In the case of superfluids there exists an intrinsic 
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contribution to F, the mutual friction force Fmf, which arises from the in- 
teraction of vortex Unes with normal excitations. This force, acting on a 
single vortex line element, ia^ 

Fmf = -apsK, X [kX (Vn - Vgl)] - a' PsK X (Vn - Vgl), (4) 

where a{T,P) and a'{T,P) are the dissipative and reactive mutual friction 
coefficients, Vsi is the local superfluid velocity at the position of the vortex 
line, K is directed along the vortex line and k is the corresponding unit 
vector. In general the averaging of the non- linear force in Eq. Q) over 
volumes containing many vortex lines cannot be done without knowledge of 
the vortex configuration. Here for simplicity we assume a locally polarized 
vortex tangle,'^ where all vortex lines in the averaged volume have the same 
direction k. In this case the transformation from F^j to the coarse-grained 
force in equation (jH)) is done by replacing k with u and Vsi with Vg. Eq. ^ 
thus becomes 



+ V(/u + -^j =VsXa; + a'u x (v^ - Vn) + au x [u x (vs - Vn)] . (5) 






This equation has traditionally been used to examine small deviations of 
vortex lines from equilibrium configurations. We use it to highlight the role 
of mutual friction in superfluid turbulence. 

The viscosity of the normal component in superfluid ^He is large: v ~ 
1 cm^/s. As a result, the normal component is in practice always in well- 
defined externally imposed laminar motion. The velocity, which is induced 
in the normal component by the mutual friction force, can be ignored in flows 
with large Rcs'- Prom Eq. (^, by substituting F with the mutual friction 
force, we find that Vn/vs ~ {Hi/v)Res (in the intermediate temperature range 
where p^ ~ Ps and a ~ 1). Thus the influence of the mutual friction force 
on the normal component can be ignored up to relatively fast motion of the 
superfluid component with Re^ ~ 2 • 10^. Even beyond that velocity (which 
has not been experimentally reached yet) the motion of the normal compo- 
nent remains laminar as Re^ ~ {K/v)'^Res. (Note that these conclusions do 
not apply to superfluid ''He where u r^ k.) 

In rotating experiments Vn = fi x r and we call v = Vg — v^ the 
counterflow velocity (Fig. ^). To simplify Eq. ((5J we change to the rotating 
reference frame (where v^ = and which we assume to be inertial): 



/u + -^ j = (1 - a')(vs X cj) + aa> X [cj X Vg] . (6) 



Q / ..2 

By applying dimensional analysis to this equation, we can construct a quan- 
tity, which has similar physical meaning as the Reynolds number in the case 
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of the Navier-Stokes equation: The ratio of the inertial (1 — a')(vs x cj) ~ 
(1 — a')y^/L^ and dissipative atj x [cj x Vg] ~ aV'^ / L? terms is here 
\/q = (1 — a')/a. Similar to the role of the Reynolds number in classical 
hydrodynamics we might expect that when 1/q is large the flow is turbulent 
and when 1/q is small the flow is laminar P Unlike in the classical case, here 
only the intrinsic properties of the fluid, a{T,P) and a'{T,P), determine 
the nature of the flow, since 1/q does not depend on the externally imposed 
geometry or the velocity of the flow. 



3. PRINCIPLE OF MEASUREMENT 

None of the established methods to generate turbulence in ^He-II are ap- 
plicable in ^He-B in the intermediate temperature regime. A new approach 
is therefore essential. Fortunately compared to ^He-II, other possibilities are 
available. One major difference is that vortex formation can be experimen- 
tally controlled in ^He-B to much better degree. If remanent vortices and 
other extrinsic sources of vortex formation can be eliminated, then a high 
energy barrier prevents vortex formation via intrinsic mechanisms up to rel- 
atively high counterflow velocities, typically of order 1 cm/s. Technically the 
simplest means of generating high vortex-free flow velocities is to rotate the 
sample container. Vortex loops can then be injected in the vortex-free flow 
by means of a few different techniques and their dynamic evolution can be 
monitored with NMR. Fig. ^ classifies schematically the different rotating 
states, from the initial vortex-free state to the equilibrium vortex state. 

The most versatile injection method at present time is the Kelvin- 
Helmholtz instability of the AB phase boundary.!^ Here the injection can 
be performed at variable values of rotation Q and temperature T. This re- 
quires a two-phase sample, like that shown in Fig. |21on the left, where the 
center section of the long sample cylinder is maintained with a magnetic 
barrier field in the A phase. 

To study the transition from regular to turbulent vortex dynamics as a 
function of temperature and flow velocity, we need a three-step process: 1) 
to set up the vortex-free rotating state, 2) to trigger the Kelvin-Helmholtz 
instability, and 3) to distinguish the type of final state. For controlling the 
experiment with the NMR measurement it needs to be able 1) to verify 
that we are correctly in the vortex-free initial state, 2) to indicate that the 
injection is properly executed, and 3) to determine the number of rectilinear 
vortex lines in the final state. Surprisingly it turns out that the transition 
to turbulence is so sharp that the final state contains either a few rectilinear 
vortex lines (Fig. ^ right), if the temperature is above the transition, or 
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else the equilibrium number of vortex lines (Fig. ^ middle) . As seen from 
the NMR absorption spectra in Fig. |2on the right, the line shapes of these 
two final states differ so profoundly that the classification in regular and 
turbulent events can be performed by inspection only. This is the case if the 
rotation velocity fi > 0.5rad/s and the so-called "counterflow absorption 
peak" (CF peak), generated by the vortex-free flow, is well developed. 

In this report we use the KH instability as injection method. It is trig- 
gered with a small step increase in rotation Afi which is typically a few 
percent of the total rotation. This means that in our measurement the ex- 
ternally controlled conditions are maintained constant, since practically even 
the rotation velocity can be regarded as constant throughout the measure- 
ment. This amounts to an important simplification of the analysis. The 
time-dependent evolution of the injected vorticity lasts from tens of seconds 
to minutes, depending on temperature and rotation velocity. During this 
transient period NMR signals can be measured which yield information on 
the dynamics. One such measurement is the determination of the flight time 
Tp of the vorticity from the AB interface to a detector coil,'^ as illustrated 
for the lower half of the sample container in Fig. [2j Another measurement in- 
volves the configuration in which the vorticity propagates along the rotating 
column .1=1 

The injection of vortex loops can be accomplished also by other means. 
In Ref. El neutron irradiation was discussed for injection. In Refs. and 
1191 a sudden burst of vortex loops through the orifice on the bottom of the 
sample space was observed to start the turbulence. 

3.1. Sample Setup 

The liquid ^He sample and the NMR measurement have been sum- 
marized in Fig. 12 Additional details have been explained in Refs. ISU] and 
1121 The long sample tube is surrounded by three superconducting end- 
compensated coil systems. These supply the two constant field regions at 
both ends of the tube for NMR polarization and the barrier field to maintain 
the A phase in the center. If the barrier field H]^ exceeds the thermody- 
namic B^A transition field Ha^b, two AB interfaces are formed symmet- 
rically around the center of the sample tube. Their shapes and locations 
depend on the current /b in the barrier magnet and on the externally con- 
trolled variables T, P, and Q^^ As an example, two locations of the AB 
interfaces are shown in Fig. [21 For the measurement of the turbulent tran- 
sition the locations of the AB interfaces are not of importance. However, to 
determine the flight time rp the distance d has to be accurately known. It is 
obtained from the calculated field profile of the solenoidal magnet system. 

The continuous wave NMR measurement is performed at constant fre- 
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Fig. 2. ^He sample and NMR measurement. (Left) The sample container is 
a quartz tube with diameter 6 mm and length 110 mm. It is separated from 
the rest of the liquid '^He volume with a partition disc at the bottom of the 
sample space. In the disc an orifice of 0.75 mm diameter provides the ther- 
mal contact to a liquid column which connects to the sintered heat exchanger 
on the nuclear cooling stage. Three superconducting coil systems produce 
independently axially oriented magnetic fields along different sections of the 
tube: at both ends are the polarization fields for NMR detection and in the 
middle the barrier field for stabilizing '^He-A. The dashed curves, approxi- 
mately symmetric about the middle of the sample tube, show the locations 
of the AB interfaces in the gradient of the barrier solenoid when its current is 
lb = 4 A (short section of A phase with parabolic contour of AB interfaces) 
and lb = 8 A (longer section of A phase with almost flat interfaces).'^ (Right) 
Two NMR absorption spectra in rotation. In the equilibrium vortex state 
the line shape borders steeply on the right to the Larmor field at 21.0 mT, 
while at low values of field sweep it extends to the temperature dependent 
cut off of the flare-out texture, here at 19.7 mT. In the vortex-free state 
the line shape displays a large "counterflow peak". Its maximum is located 
at the cut-off value of the "flare-out texture" . Both line shapes have been 
measured at the same temperature and have equal integrated absorptions. 
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quency /, using a linear sweep of the axially oriented polarization field H. 
The two absorption spectra in Fig. |21 display the extremes of different pos- 
sible line shapes. In fact, to monitor the transition to turbulence, they turn 
out to be the only two line shapes needed. The large peak in the vortex-free 
state is produced by the orienting effect of the counterflow on the order pa- 
rameter texture and is called the counterflow peak. Its shift from the Larmor 
field (at 21.02 niT) is used for temperature measurement and its height is a 
sensitive function of the number of vortex lines. 

When a cluster of rectilinear vortex lines is formed, the CF peak height 
is reduced. The reduction is linear with vortex number A^ in the limit N <C 
iVeq, if all other variables remain constant. At large values of N the reduction 
becomes nonlinear and finally, if N is continuously increased until it reaches 
the maximum possible value « N^q, the spectrum looks entirely different: 
Here much of the absorption is shifted to high fields and borders prominently 
to the Larmor edge at 21.0 mT. With increasing N the intensity in the CF 
peak is thus shifted closer to the Larmor edge. Therefore one may monitor 
the increase in N by recording either the CF peak height or the absorption 
at the Larmor edge. In the latter region the absorption is less affected by 
temperature changes and often it is more practical to monitor the evolution 
of the vortex expansion process at fixed polarization field value close to the 
Larmor edge. 

An example of such a measurement is illustrated in Fig. |31 at constant 
temperature above the transition to turbulence. Here the rotation fi is in- 
creased at constant slow rate across the region where the Kelvin-Helmholtz 
instability of the AB interface develops. The O increase is performed suffi- 
ciently slowly so that no dynamic effects are visible. The first critical event 
we denote with Qc- During further increase of the rotation drive the instabil- 
ity occurs recurrently every time when the counterflow at the outer sample 
perimeter reaches the critical value Vc = 0,cReS, where r = R^s is the ra- 
dial value at which the injected vorticity breaks through the AB interface 
fSec. I3.3"|l . Simultaneously a central cluster of rectilinear vortex lines starts 
to develop (Fig. fright). 

The parameters controlling the instability are roughly identical in the 
upper and lower halves of the container. The first instability event (with 
vortex-free flow in the two B-phase sections) is registered at almost exactly 
the same value of il in both detector coils, as seen in Fig. |21 Nevertheless, 
the instabilities of the two AB interfaces in Fig. 21 are not coupled, but 
occur independently. Thus the subsequent instabilities in Fig. El are not 
synchronous, since the number of vortices injected in the B-phase sections 
at each event varies randomly. Also in the temperature regime where the 
transition from regular to turbulent vortex expansion takes place, it is quite 
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Fig. 3. Absorption response during slow linear increase of rotation while a 
sequence of the ~ 10 first Kelvin-Helmholtz instability events is traversed. 
(Top) The measurement is performed at high temperatures where vortex 
motion is regular. Here the number of vortex loops injected in the B phase 
equals the number of new rectilinear vortex lines which are added to the B 
phase sections per instability event. The measurement monitors the increase 
in rectilinear vortex lines per event. The increase is proportional to the 
height of the discontinuous steps in the linear staircase patterns and, as seen 
from the step heights, varies randomly per event. The NMR absorption is 
measured at constant polarization field at a value just below the Larmor edge 
{e.g. for the bottom spectrometer H = 20.97 mT , see Fig. |2I). (Bottom) 
The rate of rotational acceleration is 2 • 10~^rad/s^. 



possible to have a regular event in one half and a turbulent event in the 
second half. As the two ^He-B samples are independent, one can have, for 
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instance, vortex-free flow in the upper B-phase section {N = 0) and roughly 
the equilibrium number of vortices in both the A-phase and the lower B- 
phase sections {N = Neq). 

As Qc is almost equal for both halves of the container, the first critical 
event can be triggered with a small step increase in fi at both AB interfaces 
simultaneously. This triggers the injection of vortex loops into both ^He-B 
samples and their dynamic evolution towards the respective final states is 
then recorded independently. This is our normal measuring procedure when 
searching for the transition between regular and turbulent vortex expansion. 
The top and bottom halves of the container give thus similar information, 
one half would be sufficient to identify the transition, but two halves gives 
better statistics. For the initial diagnostics, to understand the critical char- 
acteristics of the different sources of vortex formation in a new sample tube, 
the two independent NMR detection systems are most useful. Four different 
quartz tubes were examined in this setup. Their properties turned out to 
vary substantially with respect to vortex formation in the absence of the bar- 
rier field and the AB interfaces. This means that at present the cleanliness 
and the properties of the fused quartz surfaces are not in reliable control. 
The fourth and last sample tube was by far the best, with the highest overall 
critical velocity. Most of the results of this report were measured with this 
tube. 

3.2. NMR Measurement 

The classification of the vortex expansion events is based on a mea- 
surement of the number of rectilinear vortex lines in the final state. The 
vortex line number N can be accurately determined, if needed, from the 
linear reduction in the CF peak height as a function of N. It turns out that 
precise measurements of N are not needed. Because of the narrow width 
of the transition regime between regular and turbulent vortex dynamics at 
high rotation velocities, we need to distinguish only between two states, one 
with a small vortex cluster (N <^ -^eq) at temperatures above the transition 
and one with the equilibrium vortex state (A^ = Agq) at temperatures below 
the transition. In Fig. [21 the two line shapes with very different appearance 
exemplify these two cases. 

The large difference between the two line shapes in Fig. IJJis caused by 
the textural orienting influence of the macroscopic counterflow. In rotating 
experiments the textural interaction can be conveniently manipulated with 
the rotation velocity 17. The velocity of the counterflow v{N, fi) is given by 

v{N, n)=Vn-Vs = {n- n^)r = nr- K,N/{2-Kr) (7) 

within the vortex-free annulus, which surrounds the central vortex cluster in 
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Fig. 4. Spectra of ^He-B during continuous cool down at constant rotation. 
(Main panel) With decreasing temperature the CF peak is shifted further 
from the Larmor value (at 21.0 mT) and its height is reduced, since the 
total absorption is proportional to the static susceptibility xb (T, P, H) . Be- 
low 0.4 Tc the textural orienting interaction from the CF vanishes rapidly 
with decreasing temperature, as shown by this series of CF line shapes 
recorded with the bottom spectrometer at constant Q. and N = ttB? 20.^/ k < 
40 vortex lines. All spectra extend with finite (but small) absorption up to 
the Larmor value, but the line shapes above the three bottom examples have 
been truncated. (Inserts) Spectra from the top and bottom spectrometers 
measured simultaneously at the lowest temperature. In both inserts three 
different situations have been compared: a) sample not rotating (fi = 0), b) 
sample filled with the equilibrium number of vortex lines, and c) continuous 
cool-down in rotation at constant Vl with a small stable vortex cluster. 
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Fig. n (right), i.e. in the region Ry^i^^/0, < r < R. Here the rotation veloc- 
ity Qy{N) is an experimentally useful quantity which denotes the rotation 
at which the cluster of A^ lines fills the whole cross section of the sample, 
ie. the rotation velocity at which a given number of lines N is in the equi- 
librium state. No changes in the vortex-line number occur if one sweeps O 
between the limits Ov(A^) < $7 < Qc, where il.c(T,P) is the relevant critical 
velocity, at which more vortices start to form. By changing Q, one scans the 
magnitude of the textural interaction ternP21 

i^hv = -^5/9s(i-v)2, (8) 

where the superfluid density anisotropy 6ps = Psi. — Ps\\- The applied mag- 
netic field induces an orbital momentum in the direction li = Raiha where 
the unit vector h = H/H specifies the direction of the applied field H, and 
R{n, 9) is a rotation matrix in the B-phase order parameter (which minimizes 
the dipolar spin-orbit interaction in the equilibrium state). Since 5ps > 0, to 
minimize the free energy F^^ in the equilibrium configuration of the texture, 
the azimuthal CF velocity v attempts to orient the orbital anisotropy axis 
1 parallel to itself. Since the magnetic field induced distortion is of axial 
symmetry, in first order the anisotropy 6ps oc H^. In addition F^v oc Q^ (as- 
suming A^ ~ 0) and thus with increasing rotation the orienting effect from 
the CF becomes comparable sequentially with other textural interactions. 
In the axially symmetric case (when the magnetic field is oriented along the 
symmetry axis of the sample cylinder, h || z || fi) this leads to a sequence 
of different axially symmetric textures, which are generally known to be of 
"flare-out" configuration;^- 

The density anisotropy 6ps{T,P) is only weakly temperature depen- 
dent, when T > 0.4 Tc. This is the regime where all NMR measurements on 
vortices in rotating '^He-B have so far been performed and where the measur- 
ing techniques are well established. Below 0.4 Tc however, 5ps drops rapidly 
with decreasing temperature.-'^ The evolution of the NMR absorption line 
shape during cooling (in otherwise constant conditions) is illustrated in the 
main panel of Fig. ^ With decreasing temperature the peak heights rapidly 
drop, the line shapes undergo changes owing to textural transitions, stand- 
ing spin-wave resonances grow in amplitude and decorate the line shape, 
and finally the orienting influence of the CF on the texture vanishes, the 
spin wave amplitudes die away, and the NMR shifts saturate below 0.35 Tc 
(bottom spectrum) . The textural transitions as a function of the magnitude 
of the CF interaction have been studied at temperatures above 0.4 Tc 
in Ref. 1241 Broadly speaking the transitions in Fig. |^ resemble those which 
occur at higher temperatures while is reduced at constant temperature. 
An example is the second CF peak at maximum possible NMR shift which 
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starts to grow in amplitude at T < 0.42 Tc. At still lower temperatures this 
second peak breaks into a twin peak structure while the main CF peak at 
smaller NMR shift is already declining in amplitude. The doubling of the 
peak at maximum possible NMR shift has not been observed before and no 
explanation exists for it. Finally below 0.35 Tc the effect of the counterflow 
is not visible in the line shape and the final example on the bottom of Fig. |^ 
is similar in appearance to the equilibrium vortex state spectrum in Fig. [2j 

The two inserts to Fig. |3 demonstrate that at the lowest temperature 
the NMR spectrum is completely insensitive to rotation: The line shapes 
are identical irrespective whether the sample (a) is not rotating at all (O = 
0), (b) rotating in the equilibrium vortex state (A^ = Agq); (c) or with 
vortex-free flow around a small vortex cluster (f^v = 0.050 rad/s). The NMR 
polarization field of the top spectrometer is 40 % higher than that of the 
bottom spectrometer. In the upper insert of Fig. H] the textural magnetic 
healing lengthP^ is of order ^h ~ 3 -R and the line shape is that of a regular 
flare-out texture. In the lower insert at lower NMR field ^h oh 1/H is larger 
and the line shape is more distorted. The spectra from the top and bottom 
spectrometers in these two inserts have been measured simultaneously in 
zero barrier field (i.e. in the absence of an A-phase layer) at the lowest 
temperature to which the sample cools. This temperature is limited by the 
residual heat leak into the ^He-B sample and its thermal contact to the 
refrigerator. Owing to the saturation of the NMR shifts, and to a lesser 
degree the uncertainty in the value of the Larmor edge (i.e. the line width 
from the residual magnetic field inhomogeneity of the NMR polarization 
magnets, FWHH r^ 0.8 kHz), with the presently available techniques the 
lowest temperature can only be assigned an upper limit T < 0.35 Tc. 

An interesting further feature in the line shapes of Fig. 0] are spin- wave 
resonances in the two signals just above that on the very bottom. In contrast 
to the local oscillator model, which yields the overall line shape of the flare- 
out textures, spin-wave absorption modes represent coherent spin precession 
over global regions with a smoothly varying texture.'^ The magnetic healing 
length of the texture ^h grows with decreasing temperature and becomes 
finally comparable to the sample radius R. This causes gradients in tex- 
tural orientations to become more and more gradual. Thus the spin-wave 
resonances grow in prominence with decreasing temperature. As seen in 
Fig. m both the intensity and the width of successive spin-wave resonances 
are approximately constant across the entire absorption regime of the sig- 
nal. In fact, theory predicts that their absorption intensity grows towards 
low temperatures and finally saturates, similar to NMR shifts. However, 
in the signal corresponding to the lowest temperature in Fig. |1] spin-wave 
resonances are not visible. It is believed that here, owing to the absence of 



A. p. Finne et al. 

dissipation and rigidity of the order parameter texture, spin-wave resonances 
are not displayed by our measurement. 

The most convenient thermometer in these measurements is the hori- 
zontal shift of the CF peak from the Larmor value. The shift is calibrated 
in a slow continuous warmup, which is repeated at varying warming rates, 
against two external thermometers. These are rigidly fixed on the nuclear 
cooling stage: a capacitive ^He melting pressure gauge and a pulsed NMR 
thermometer with a platinum wire brush as sample. Our temperature cal- 
ibration of the NMR shift is generally in agreement with those reported in 
Refs. I2SI and EHl However, below 0.7 mK the external thermometers lose 
contact with the ^He-B sample and a smooth extrapolation of the NMR 
shift has to be used as temperature scale. Finally, as seen in Fig. |1] below 
~ 0.35 Tc, the NMR shifts saturate and become useless for thermometry. 

To summarize, in the present type of NMR measurements the experi- 
mentally vital information is obtained via the image which the NMR pro- 
vides about the order parameter texture, primarily through the change in 
the magnitude of the counterflow interaction in Eq. (jHJ as a function of 
the number of rectilinear vortex lines N at constant Q. These techniques 
work well at temperatures above 0.4 Tc, so that single-vortex resolution is 
achieved above 0.7 Tc. However below 0.35 Tc the texture-based NMR tech- 
niques break down. This is most unfortunate for further efforts to study the 
properties of turbulence in the zero temperature limit. Clearly T > 0.35 Tc 
is not yet the regime of ballistic quasiparticle motion and vanishingly small 
mutual friction. At present it is not clear if other NMR techniques can be 
developed for ^He-B to measure the number of rectilinear vortex lines or the 
vortex density at the lowest temperatures. 

3.3. Vortex Loop Injection 

The classic Kelvin-Helmholtz (KH) instability is a generic phenomenon 
of fluid dynamics which occurs at the interface between two inviscid fluid 
layers, which move with respect to each other in a state of relative shear flow. 
If the difference in tangential velocities is small, the interface remains calm, 
but at some critical value interfacial waves are formed. The AB interface 
instability is its first superfluid example.^ 

^He-A and B behave very differently in rotation: The A phase has a 
low critical velocity and remains close to the equilibrium vortex state at all 
times, while the B phase is vortex-free before any instabilities have occurred. 
To maintain this situation, a new structure of vorticity is formed on the A- 
phase side of the AB interface, a surface vortex layer. It is confined to the 
interface and gives rise to the tangential velocity difference across the phase 
boundary.l^ZI The instability manifests itself as a sudden event in which some 
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Fig. 5. A sequence of KH instability events with increasing fi. The ab- 
sorption amphtude close to the Larmor value has been recorded with the 
bottom spectrometer. The sloping dashed line is a fit through the insta- 
bility points and is proportional to A^ = 7ri?^g2(r2 — i^c)/n- (Bottom) The 
discontinuous CF velocity as a function of Q at r = R^ff-. v =| UsB — '^n | = 
i^ReS — nN/{2'KRes)- The horizontal dashed line is equivalent to the sloping 
dashed line in the upper panel and defines the mean critical velocity Vc- 



vortices from the surface layer break through the interface and start to evolve 
in the vortex-free B-phase counterflow. 

The critical CF velocities across the interface follow the general formula^Sl 



2 /OsA(WsA - V^f + - Psb{VsB - Vnf = \/(^AbF . 



(9) 



A. p. Finne et al. 



I I I I I I I I 
r= 0.774 Tc 
P = 29.0 bar 
/b = 4.0 A 




5 10 15 20 25 

Number of B -phase vortex lines/KH event 

Fig. 6. Histogram of 214 measured shear-flow instability events of the AB 
interface, plotted as a function of the number of rectilinear B-phase vor- 
tices AA^ formed per event. The smooth curve is a guide for the eye of a 
differential probability distribution. 



where the magnetic restoring force on the AB interface is F = -^ (xA — 
XB)V(-ff^) at Hy^ = Hab- Here xa(-P) and xniT, P,H) are the suscep- 
tibilities of the A and B phases, axB{T,P) the surface tension of the AB 
interface, and Hab{T,P) the critical magnetic field of the AB transition. In 
our experiments v^a ~ ^n and v^b — Vn = ^cRcS, where Q,c is the first critical 
rotation velocity of the KH instability and R^s the effective radial value for 
the instability. The wave vector of the mode of interface waves, which causes 
the instability, is A; = \J Fja . Measurements in widely diff'erent conditions 
agree with Eq. @ without adjusted or fitted parameters W^ 

A repetition of Fig. |21 at higher signal-to-noise resolution is shown in 
Fig. [5] for the bottom section of the sample. It illustrates the experimental 
consequences from Eq. @ , the signatures of the KH instability. A sequence 
of instability events has been recorded here in slowly increasing rotation at 
high temperatures where the number of vortex loops injected in the B phase 
equals the number of rectilinear vortex lines which are formed after each 
instability event. The signal monitored in the top panel is the absorption 
close to the Larmor edge where the step-like increase measures the number of 
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new rectilinear vortex lines which are added to the lower B phase section, i.e. 
the vertical discontinuous step height is proportional to the increase in vortex 
number. The simplest technique to assign the correspondence between step 
height and vortex number AA^ is to perform a least squares fit with the 
constraint that the values for AA^ have to be integral numbers. The order of 
magnitude for the step height signal is 2 ^V/vortex line. This figure refers to 
the voltage measured at the output of our cryogenic MESFET preamplifier 
(with a gain ~ 10) which follows the superconducting tank circuit. With 
this fit we get the vortex line numbers which are denoted in Fig. [5] next to 
each step. 

The location of the instability events on the axis has to fulfill the 
requirement that the critical KH velocity is constant, i.e. the corner points in 
the staircase pattern fall on one line. This gives the dashed linear fit with the 
slope 1.27//V/(rad/s). Combining the two numbers so far, we get AN/Ai} = 
671 lines/ (rad/s). This means that when AA'^ new rectilinear lines are added 
to the central vortex cluster, the counterflow velocity along the effective 
perimeter decreases by AQ ■ R^s = kAN/{2ttR(.q). From this requirement 
we obtain the effective radial value i?eff ~ 2.6 mm for the instability, which 
should be compared with the sample radius R = 3.0 mm. 

A lot of similar examples confirm the above value for the effective radius 
i?cff < R- Also more elaborate measurements, which measure the number 
of vortex lines formed in each instability event with a calibrating procedure, 
give the same value of i?eflf in the above regime of O, T, P. A result from 
this latter type of precision measurements was reported in Ref. |27| What 
then is the explanation for i?cff? The vortices which create the staircase 
pattern in Fig. [51 are detected by the pick-up coil in the bottom of the 
sample tube 29 mm below the AB interface. The flight time for the vortices 
from the interface to the top edge of the coil is 5 s. At high temperatures the 
expansion of the vortices is deterministic. The fixed flight time introduces a 
5 s delay before the information about the events at the interface starts to 
build up in the pick-up coil. In Fig. [3] the rotation drive is a linear ramp with 
the acceleration rate dil/dt = 2.5-10"^ rad/s^. This means that every 6 s the 
counterflow velocity increases at the effective perimeter by the equivalent of 
one circulation quantum. Comparing these two numbers, we interpret the 
staircase pattern to arise from the events at the interface and not to be 
blurred by the motion of the vortices along the rotating column. In this 
case i?eff is the radial distance at which the KH instability occurs and the 
vorticity is injected across the AB interface from the A phase into the B 
phase. 

The fact that we measure a definite value for Res allows us to make 
a guess of the initial configuration of the injected vortex loops immediately 
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after injection, in the situation from where their expansion in the vortex-free 
B-phase counterflow starts. This would be a loop of length > i? — Rcs ~ 
0.4 mm whose one end sticks out of the AB interface at r = Res and the 
second ends on the cylindrical side wall at r = R in the B phase section. 

The second important question about the injection is the number of 
such loops. In Fig.lSJwe have a small sample of eight instability events, with 
the number of vortex lines AA^ formed per event marked next to each event. 
It is seen that AA^ is a small random number and on an average AA'^ ~ 7. 
In Fig. El a histogram of AA^ for a large sample of 214 events is shown. Such 
a probability distribution can only be measured at high temperatures with 
regular vortex motion. It has not been carefully measured as a function of 
the external variables T, 0, P, but there are no indications that it would 
be rapidly changing in the conditions used for vortex loop injection in this 
work. Therefore we assume that the distribution behaves in a continuous 
manner across the transition from regular to turbulent dynamics and is of 
similar shape in the turbulent regime. 

The average of the probability distribution in Fig. [SJis 11 vortex lines 
per event. This agrees with what one would expect for the number of A- 
phase vortex quanta in the surface vortex layer which cover one trough, or 
half a wave length, in the wave pattern of the interface ripplons according 
to Eq. 0. In the A-phase section, assuming solid-body rotation with zero 
critical velocity, there are A^ ~ ^^cff "^^c/f^ single-quantum vortices (k = 
h/2m^) covering the AB interface which all flare out radially perpendicular 
to the cylindrical wall. Measured along the effective perimeter there are thus 
Jlc-Reflf/^ circulation quanta per unit length and thus in one corrugation 
AA^ ss {■KVc)/{kK,), where Vc = i^cReS- From Eq. @ this is seen to be 
AA'^ ~ {2Tra)/{KVcPsB), which gives 9 vortex lines in fair agreement with 
the measured average in Fig. El (fc = 0.39 cm/s at 0.77 Tc, 29.0 bar, and 
/b = 8.0A). 

As noted in the context of Fig. 01 the first critical event, before any 
vortices have been injected, takes place almost simultaneously in both B- 
phase sections. Therefore, when the critical rotation velocity ilc is known, 
the two events can be triggered together with a small step increase Ail. = 
0.050 rad/s in the manner shown in Fig.d This procedure allows to measure 
the flight time for the vorticity to propagate from the AB interface to the 
detector coil. Moreover, it allows to set up the NMR measurement so that 
transient signals can be captured efficiently when the vorticity front travels 
through the detector coils. In contrast to Fig. |21 in Fig. |3 the evolution of 
the CF peak height after injection has been monitored at low temperatures 
where the expansion of the injected vortex loops is turbulent. 

What happens at injection? What is the configuration in which the 
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Fig. 7. Response of tlie CF peak lieigiit to a triggered KH injection event 
in tiie turbulent temperature regime. (Bottom) The instability is triggered 
by increasing fi from 1.31 to 1.36rad/s in 2 s. With vortex-free flow in both 
B-phase sections, the instability at both AB interfaces falls between these 
two limits. (Top) The trigger step increase AQ is instantaneously registered 
as an increase in the CF peak height. A rapid decay of the peak height 
starts when the front of propagating vorticity reaches the closest end of the 
detector coil and approaches zero when the front leaves the far end of the coil. 
This signal decay corresponds to the removal of the large-scale azimuthal 
counterflow just behind the propagating front. The NMR polarization fields 
are maintained here at constant value at the location of the maximum of the 
CF peak of the initial state {e.g. at 19.37 mT for the bottom spectrometer 
in Fig. E}. 
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vorticity expands along a rotating superfluid column? The answers to these 
questions are still discussed; we are not going to delve into details here. 
Clearly the first turbulent proliferation of vorticity takes place immediately 
at or after injection. If the rotating vortex-free column is infinitely long, it is 
unrealistic to expect that turbulence would fill at some point of the process 
the entire column. Rather it could be expected to expand into the vortex- 
free section in a form which after some time travels along the column in a 
time-invariant average configuration. Since the flight time of the vorticity 
grows as rp = d/{aO,R) with decreasing friction a, the propagating front 
moves slower with decreasing temperature and has more time to settle down 
in such a configuration. Thus at 0.40 Tc the transient signals do not appear 
to depend on the distance d. This conclusion is reached from monitoring the 
transient signals, which the propagating front induces in the two detector 
coils when turbulence is started in random locations along the sample with 
a neutron capture process (in zero barrier field and no A-phase) .l^ With this 
picture in mind, we take a second look at the timing in Fig. [3 

The split-half detector coils, whose axes are oriented transverse to the 
sample column in the setup of Fig. 121 have a sharp cutoff in sensitivity at the 
edge of the coil. Thus the flight time from the AB interface to the detector 
coil is bracketed between the instantaneous increase in CF peak height from 
the trigger and the arrival of the vorticity front at the closer end of the 
detector coil (from where the decay of the CF signal starts in Fig. [T]). In 
Fig. 121 this flight time of 41s corresponds to a distance d ~ 34 mm (bottom 
section of sample). The front is expected to travel through the 10mm long 
detector coil during the next 10/34x41 s = 12 s. From Fig. Owe note that the 
CF peak height of the bottom spectrometer vanishes in ~ 14 s. Comparing 
these two numbers, this means that the macroscopic azimuthally flowing 
counterflow is removed immediately behind the traveling front of vorticity, 
i.e. the vorticity behind the first front has the equilibrium polarization of 
solid body rotation in the direction z of the column. Other conclusions about 
the structure of the front in the turbulent temperature regime can be drawn 
from the transient signal measured at the Larmor edge, i.e. by monitoring 
the build up of the vorticity instead of the CF decay.'^ 

Figs. and d thus characterize different aspects of how the vorticity 
expands in the rotating column at high and low temperatures. The most 
fundamental difference is in the final states: In Fig.Oonly a small number of 
rectilinear vortex lines is created in each sequential KH event, while in Fig. d 
the first event fills the vortex-free B-phase section with the equilibrium num- 
ber of vortex lines. We have measured extensively the first critical velocity 
r^c of the KH instability as a function of T, I\,, P. These measurements give 
in the (0, T) plane continuous trajectories,!^ in which the transition to turbu- 
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lent vortex dynamics is not visible. There is no indication that the injection 
process itself would influence which final state is formed, when the tempera- 
ture is well above or below the transition. Rather, all measurements suggest 
that the KH instability is a predictable injection mechanism with continu- 
ous properties and with a reasonably narrow distribution in the number of 
injected vortex loops. 

The measurement of flight time rp in Fig. [3 gives the propagation ve- 
locity of the vorticity along the rotating column: vi^^ = d/r-p. It turns out 
to amount to wlz = aJ^-R, since this expression gives a values in good agree- 
ment with the mutual friction measurements in Ref. 1101 This is the value 
expected for the end point of a single vortex spiralling down along the ver- 
tical container wall in the vortex-free state. The measurement can thus be 
used to determine the dissipative mutual friction coefficient a(T,P). Sur- 
prisingly the single-vortex value for ulz applies for the propagation velocity 
both in the regular and turbulent regimes of propagation, i.e. the velocity of 
the front of expanding vorticity with the retreating vortex-free flow displays 
no anomaly at the transition.'^ 



4. TRANSITION TO TURBULENCE 

When measured with KH injection, the transition from regular to tur- 
bulent dynamics is found to be velocity independent and centered at 0.52 
- 0.59 Tc, depending on pressure.l^ It coincides with the condition q = 
a/(l — a') ~ 1. However, the transition is liquid pressure dependent, mov- 
ing with increasing pressure to higher temperatures on the normalized T/Tc 
scale and also to higher q. Thus the transition does not take place at exactly 
the same q value at all pressures and the situation is more complex than the 
simple theory in Sec. ED portrays. With KH injection the transition has a 
narrow width of 0.06 Tc at the three measured pressures of 10.2, 29.0, and 
34 bar. In this section we describe results from the transition region which 
shed more light on the KH injection mechanism. 

A close up of the transition regime between regular and turbulent vortex 
motion as a function of temperature is shown in Fig. |S1 Here the number 
of rectilinear vortex lines has been measured and is expressed as a fraction 
of the equilibrium state, i^^/Q, where Q is the critical rotation velocity 
Jlc of KH injection. The conspicuous feature is that intermediate values 
f^v/^ ~ 0.1 - 0.7 are prominently absent: the final state after injection 
includes either very few vortices or is close to equilibrium. In fact, the only 
intermediate values which have been measured are a few in the transition 
regime. This means that turbulence predominantly either switches on and 
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Fig. 8. Number of rectilinear vortex lines formed in a KH injection event 
as a function of temperature in the transition regime. The measured re- 
sult on the vertical scale gives the fraction Oy/fi, the number of vortices 
present compared to the equilibrium number. This ratio is obtained using 
a calibration procedure similar to that explained in Ref. 1121 In most cases 
the result pertains to the first injection event with vortex-free flow in the 
B-phase section, but some cases have been included where the first or even 
the second sequential event gives only a small number (as in Fig. ^ , until a 
later injection event produces a large number of vortices. The measurement 
demonstrates that the temperature width of the transition regime is narrow. 
Also the number of rectilinear vortex lines produced per injection event is 
either very small and the process is of the form shown in Fig. IHl or the final 
vortex number is close to equilibrium and the process is as in Fig. [7| In 
contrast, intermediate vortex clusters are rare. 



fills the sample with close to the equilibrium number of vortices or it does not 
and the final state contains only the vortices injected by the KH instability. 
Therefore the transition regime is not characterized during cooling by a 
gradual increase in the value of Q^/O,, but by the fact that a growing fraction 
of the expansion processes after injection become turbulent and send the 
sample in the equilibrium vortex state. 

The absence of intermediate values for the fraction Q^/Q has an inter- 
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Fig. 9. Sequence of two turbulent transitions during a slow continuous ac- 
celeration of the rotation drive. The absorption response of the bottom 
spectrometer has been followed at constant polarization field just below the 
Larmor edge. 



esting explanation. It suggests that turbulence either switches on or does not 
depending on initial conditions, as was observed to be the case in neutron 
capture induced turbulence in Ref. |^ The initial conditions are dictated 
by the details of the injection process. Although the KH instability appears 
to be a reproducible injection technique, it contains a stochastic component, 
namely the number of injected vortex loops AA^ in Fig. El and their exact 
configuration. This would mean that it is not the expansion process, but the 
injection situation, which is most influential in starting turbulence: in the 
same conditions the expansion process might proceed towards only a few vor- 
tices or to the equilibrium number and it is the stochastic variability in KH 
injection which produces the finite width of the transition. These considera- 
tions are corroborated by ongoing numerical simulations of the appropriate 
experimental situation.!^ 

A further example of the transition features is illustrated in Fig. |^ 
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It shows the NMR absorption response, measured at constant polarization 
field close to the Larmor edge, during a continuous slow increase of O. The 
beginning of the trace (0 < 1.4rad/s) shows how the absorption is shifted 
from the Larmor region into the growing CF peak (compare to the spectra 
in Fig. |2)). Here it can be seen that the CF peak starts to develop only 
when Q > 0.2rad/s. Above 1 rad/s the textural interaction from the CF has 
become the dominant contribution and from then on changes in the Larmor 
region are minor until the KH instability occurs at 1.45 rad/s. The sharp 
peak here is a transient signaP from a turbulent expansion process which 
fills the sample with almost the equilibrium number of vortices. 

During further increase of fi a central cluster of vortices with 0,^ ~ 
1.4 rad/s is formed, but otherwise a similar cycle of absorption behavior is 
repeated. Finally a second turbulent event takes place at exactly twice higher 
rotation velocity, as appropriate for the critical velocity of the KH instability 
according to Eq. Q. This latter turbulent event occurs in the annular 
vortex- free CF region which surrounds the central cluster in Fig. ^ (right). 
This example shows that (i) the central cluster formed from rectilinear vortex 
lines is stable in increasing rotation and (ii), since the second turbulent 
event occurs at the expected Q value of the KH instability, the turbulent 
event evolves outside the cluster at large radii and CF velocities, so that the 
cluster is not involved in the process. 

In Fig. ^Jresults are displayed on the final state after KH injection, as a 
function of CF velocity at constant temperature. The two plots examine the 
velocity dependence of the transition and here exactly in the temperature 
regime where the transition from regular to turbulent dynamics takes place. 
The data have been classified as (a) a regular expansion event if the final 
number of rectilinear vortex lines is in the range of the distribution in Fig. El 
as (b) turbulent if the final state contains close to the equilibrium number 
of rectilinear vortex lines, or as (c) mixed if one or a few regular events 
precede a final turbulent event. The data have been collected by moving 
as a function of the current I^, in the barrier magnet along an isotherm in 
the (0,T) plane. Since T and P are constant in this measurement, also 
the magnetic field at the interface is constant: H = H^b- On the left at 
low currents the curves start when the maximum field in the center of the 
barrier solenoid at z = reaches H^-q. From here with increasing lb the AB 
interfaces move apart along the solenoidal characteristic H{I\i,z) = H\b 
and the gradient VH \ h=h ab increases. As a result, the KH critical velocity 
also moves to higher values, which is seen both from Eq. Q and in Fig. 1101 

Among the many different measuring runs which have been performed, 
the isotherms in the two panels of Fig. ^J happen to be almost exactly in 
the center of the transition regime at 10.2 and 29.0 bar pressures. They were 
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Fig. 10. Final states after KH injection in the temperature region of tlie 
transition from regular to turbulent vortex expansion, classified in (a) reg- 
ular expansion events (open triangles), (b) turbulent events (filled squares), 
and (c) cases of one or a few regular events preceding a final turbulent event 
(triangle over a square). The curved solid lines represent Eq. @ without 
adjustable parameters .1221 The insert in the bottom left panel is a cartoon of 
the interfacial wave on the AB phase boundary. The black half circles rep- 
resent the A-phase vortex quanta covering the interface and flaring radially 
outwards {i.e. out of the plane of the paper). 
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selected to illustrate the behavior in the transition regime. The results show 
no clear velocity dependence: The mix of different symbols in both panels 
appears random. This is also the general conclusion reached in Eef. 1131 when 
all data from the KH measurements is reviewed as a whole in the $7 range 
0.8- 1.6rad/s. 

The solid curves in the panels of Fig. ^J illustrate Eq. Q and its pre- 
dictions. The panels in the middle show the behavior of the half length of 
the interface wave at the instability: ^ X = n y^OAB/F- Inserting the critical 
KH velocity ^ psB ^c — Vo'AB F, we obtain 

1a = 2.^^-1 (10) 

and therefore the trough with the A-phase vortices, which break through 
the interface, gets rapidly smaller with increasing rotation. The inter- vortex 
distance a in the surface vortex layer on the AB interface is determined by 
the tangential velocity difference: a{vsB — ''^sa) ~ n (see insert in panel on 
bottom left). The number of vortices in the trough can then be expressed 
as 

A AT -^ o '^AB 1 „ f^AB ...x 

AiV r^ — = 2tt -—- = ZTT — ^5 a (11) 

2a kPsbRcS "c k PsB 

and thus it also decreases with increasing rotation, as shown in the outermost 
panels on the right. 

So far the probability distribution in Fig. IHlhas not been measured as a 
function of $7, to check the dependence AA^ oc Q,^^. However, it is instructive 
to consider the consequences from Eq. (|11() . Our measurements at 29.0 bar 
and lb = 8.0 A as a function of temperature are a good example. First 
i^ciT) was measured using a slowly accelerating rotation drive, as in Figs. IHl 
and|5J ^c{T) turned out to increase with decreasing temperature, curving 
smoothly towards saturation at ilc ~ 1.6rad/s below 0.6 Tc. This agrees 
with Eq. @ and reflects the fact that at 8 A the barrier field in the center 
of the solenoid exceeds Hab (T ^^ 0)- From Eq. (fTT|) we see that if anything, 
then AA^ decreases towards low temperatures. The measurements were next 
repeated using the triggering method of Fig. [3 with AQ = 0.050 rad/s. Such 
a measurement proceeded as seen in Fig. |7|and resulted below the transition 
regime < 0.53 Tc always in a turbulent event. In view of Eq. (|1H) this seems 
to mean that even with fewer injected vortex loops the first KH instability 
leads invariably to turbulence, or that even a small number of injected loops, 
such as AA^ < 3 (Fig. EJ, has to suffice to start turbulence. This we take as 
evidence for the importance of the Kelvin- wave instabilitjff^of a single vortex 
as the initial source for the proliferation of vorticity and as the mechanism 
which here starts the turbulence. 
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Finally we note in passing two features which are related to any future 
extension of these measurements: (i) With the present barrier magnet the 
trajectories of the KH instability do not extend much below or above the 
velocity interval shown in Fig. 1101 For lower velocities a barrier magnet 
with a weaker gradient is needed, for higher velocities it should be steeper 
(Eq. ©)• (ii) We have performed measurements with KH injection also 
at zero pressure, to span the entire range in particle densities. However, 
at low temperatures below ~ 0.5 Tc a new interfering feature appears, the 
continuous proliferation of vorticity when curved vortex lines are present and 
undergo the Kelvin-wave instability.'^ It then becomes difficult to achieve 
vortex-free rotation to sufficiently high velocities, as required for the KH 
instability. This difficulty grows in prominence at low pressures, especially 
at zero pressure. The Kelvin-wave instability starts to fill the sample at a 
seemingly constant rate with rectilinear vortex lines already at low rotation 
and proceeds until close to the equilibrium vortex state. The source of the 
problem are existing curved vortex lines, which at = are remanent 
filaments, left over from earlier rotations. These are difficult to avoid at low 
temperatures, since the last one or two remanent vortices take exceedingly 
long times to annihilate in our long sample tube in standstill (Q = 0). 
Rotation in the reverse direction is not helpful since the required velocities 
are of order ~ 10^^ rad/s and even then more often than not lead to a 
reorientation of the filament (a reversal in the orientation of k.) and thereby 
via the Kelvin- wave instability to more vorticity. The long life time of vortex 
filaments at the lowest temperatures in zero applied flow we noted when the 
example of a single precessing vortex from Ref. I^was described in Sec. IlI 



5. CONCLUSIONS 

We have described the techniques for investigating turbulence in ^He- 
B in the temperature range T > 0.4 Tc. Both the methods to generate 
turbulence and to record it non-invasively are different from what can be 
done in He-H. In the temperature region of these measurements the vortex 
mutual friction of ^He-B changes rapidly and proves to control the character 
of the vortex dynamics. A sharp transition from regular vortex number 
conserving dynamics at high temperatures to turbulent dynamics at low 
temperatures is observed at T < 0.6 Tc. The best injection mechanism of 
vortex loops into vortex-free counterflow of '^He-B has been found to be the 
Kelvin-Helmholtz instability of the magnetically stabilized interface between 
^He-A and ^He-B. In this process only the number of injected vortex loops is 
not in good control, which is believed to cause the observed finite transition 
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width. 

It is to be expected that further work on turbulence in ^He-B wih al- 
low to place turbulence in ^He-II, the only previously known example of 
superfluid turbulence, into a wider context. In contrast to ^He-II, where 
in most cases at higher flow velocities one has to deal with the mutual- 
friction coupled turbulence in both the superfluid and normal components, 
in ^He-B experiments the normal component is practically always immobile. 
Thus a simpler example of superfluid turbulence is available for studies here. 
The present NMR technique, which relies on the influence of the vortex-free 
counterflow on the order parameter texture, stops working below 0.35 Tc. 
At present it is not known if other NMR techniques can be developed to 
investigate turbulence in the zero temperature limit. 
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